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Abstract 

The ring-shaped Hartmann potential V = rja^eQ — :pT^^^ was introduced 
in quantum chemistry to describe ring-shaped molecules like benzene. In this arti- 
cle, fundamental concepts of supersymmetric quantum mechanics (SUSYQM) are 
discussed. The energy eigenvalues and (radial) eigenfunctions of the Hartmann 
potential are subsequently rederived using the techniques of SUSYQM. 
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1 Introduction 



In 1972, an exactly solvable ring-shaped potential was introduced by H. Hartmann 
The Hartmann potential is given by the following expression, 

V = 7]a tol (1) 

where 

fi^ 1 /ie^ 

ao = — ^ and = tj- (2) 

fi is the particle mass, 77 and a are positive real parameters which range from about 1 to 10 
in theoretical chemistry applications |^ and r, 9 are in spherical coordinates. Following 
the exact solution of the Schrodinger equation given by Hartmann IQ , alternative methods 
of solutions in spherical coordinates and squared parabohc coordinates 0, ^ H, ^ |^ 
had been given. In this article, the author presents another alternative method of solution 
in spherical coordinates using supersymmetry (SUSY). 

The concept of supersymmetry (SUSY) has been used in particle physics in the past 



two decades [§, |TOl- It was discovered in 1971 by Gel'fand and Likhtman [|lT|. Simply put 



supersymmetry is a symmetry which relates fermionic and bosonic degrees of freedom. 
At present, particle physicists believe that it is an essential ingredient in unifying the four 
fundamental forces in nature namely the electromagnetic, weak, strong and gravitational 
interactions. 

Supersymmetric theories of the four fundamental interactions entail the presence of 
SUSY partners which have the same mass as their corresponding ordinarily observed 
particles.^] Unfortunately, these have not been observed in nature. To make sense out 
of this experimental fact, theorists believe that SUSY must be "broken" at ordinary 



energies. The search for a mechanism to break SUSY led Witten |jT2[ in 1981 to study 
SUSY breaking in the simplest case of SUSY quantum mechanics. In fact, studies in 
SUSYQM during its early years, were confined solely for understanding SUSY breaking. 

However, it was eventually discovered that SUSYQM can have interesting applications 
besides its use in the study of SUSY breaking. At present, it has found its way in many 
areas of physics including atomic physics, statistical physics, nuclear physics, etc. p!3| . 
Through the present article, the author hopes to contribute to the utilization of SUSYQM 
in theoretical chemistry. 



^For example, the SUSY partner of the electron is called a selection while that of a photon is a photino. 



2 



One simple use of SUSYQM is in obtaining the eigenvalues and eigenfunctions of 
the Schrodinger equation [0, |T^. A two-dimensional supersymmetric solution of the 
Hartmann potential in "square" parabolic coordinates 0] had been proposed. In the 
present article, the author formulates a one-dimensional SUSY solution similar to that of 



the hydrogen atom |jT^. This avoids the complications of a two-dimensional SUSYQM 
formulation. Using SUSY, it is shown that the eigenvalues and (radial) eigenfunctions 
can be obtained. The techniques employed here are largely due to Schwabl [|T^. 

In section |, a pedagogical introduction to SUSYQM is developed. Only the concepts 
and equations which are essential for the present paper are presented. 

Section |^ gives a detailed account of how SUSYQM can be used to obtain the eigen- 
values and (radial) eigenfunctions of the Hartmann potential. The discussion is heavily 
supplemented by energy level diagrams to facilitate the understanding of the hierarchy of 
hamiltonians and how the SUSY operators Aj^ are used to obtain the eigenfunctions. 

Some conclusions are given in section 0. 
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2 Supersymmetric Quantum Mechanics 



As mentioned in the Introduction, SUSY was first applied to particle physics, whose 
language is quantum field theory. In quantum field theory, a particle is represented by 
a component field y?, and its dynamics is described by a lagrangian density C{ipi,d^ipi) 



where 



1^ V 



. The word "supersymmetry" was originally used to describe the 
symmetry which transforms a field Lp to another field if) whose intrinsic spin differs from if 
by ^h. In SUSYQM, which will be described here, we will use the term "supersymmetry" 
in a more general sense. It will be used to denote systems which can be described by the 
SUSY algebra in supersymmetric field theory. 

SUSYQM [|T2], |T^ is characterized by the existence of the charge operators Qj, where 
i = 1, 2, . . . , such that they obey the SUSY algebra (denoted by sqm(A^)), 



[Qi , Hg 







(3) 



where Hgs is the supersymmetric Hamiltonian, { } and [ ] are anticommutator and com- 
mutator respectively. We consider only sqm(2) with charge operators Qi and Q2 and 
construct the linear combinations 



and 

From equations ^ and |^, The SUSY algebra is then 



Q = ^{Qi + iQ2) 



'I 



with 



and 



0. 



The above SUSY algebra can be realized by letting 



Q 





A- 









(4) 



(5) 



(6) 



(7) 



where 



From equations ^ and |^, the supersymmetric hamiltonian is 



Hi 
H2 



(9) 
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where 



Hi = A+A- 



and 



A-A-^ 



(10) 



The hamiltonians Hi and H2 are said to be "supersymmetric" partners of each other. 

The hamihonian of the Schrodinger equation can always be factorized in the form of 
equation [l^. Consider the hamiltonian |]13| 



Hi 



(11) 



such that 



(1) 



+ vAx) 



^(1) = ^(1)^(1) 



where Vi{x) is chosen such that the ground state ip^^^ has energy equal to zero 



hamiltonian in equation 11 can be put in the form of equation M by letting 



1 / d 



and 



At 



\dx 

provided that the "superpotential" Wi satisfies the Ricatti equation 



(12) 
The 

(13) 



Vi{x) 



wl - 



dWi 
dx 



(14) 



As long as equation O has a solution Wi, the one-dimensional Schrodinger equation can be 



made supersymmetric by the construction given in equations |13|, |T^ and ^. The challenge 
then in using SUSYQM techniques is not in the mechanics of the construction of just 
any SUSY hamiltonian, but in finding a suitable superpotential (or Vi{x)) to construct a 
SUSY hamiltonian which will be relevant to the problem at hand. It is a common practice 
to choose or pose as an ansatz the Wi to solve a physical problem |14, l7 . 
The SUSY partner of Hi^ namely H2 is then given by 



where 



V.ix) 



dx 



(15) 



(16) 



Note that H2 is altogether a new hamiltonian. An astute reader will immediately 
realize that one can repeat the procedure in constructing H2 from Hi to construct an H^ 
from H2 such that 

H2 = AtA2 (17) 
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with 



' V2 \dx \ 



and 



At=^ i-i- + W, 



y/2\dx 



and with a new Ricatti equation 



V2{X) = - 



dW2 
dx 



W2 in equation 19 is then solved to construct A2 . H3 can then be constructed as 



^2 ^2 



where 



V.(x) 



1 



dx 



(18) 
(19) 

(20) 
(21) 



We can evidently construct a "hierarchy" of SUSY-partner hamiltonians Hi, H2, H^, . . . , 
Hn starting from Hi. 

Let us go back to the first two hamiltonians we started with namely Hi and H2. Since 
Vi (x) in equation |TI| is chosen such that its ground state wave function tp^^^^ has eigenvalue 
equal to zero, equation |l^ gives 

i/i^°i)=0 =^ A+ArV^(i)=0. (22) 

Equation ^ implies 

^r<i) = 0. (23) 

With equation |T3| and knowing Wi, equation allows one to calculate the ground state 
of Hi by solving the resulting first order differential equation. 
Consider any eigenstate of Hi, ■j/'"^^ with energy EJ\y We have 



Hl^l,^l~^ = El^r^ 



1) 



(24) 
(25) 

(26) 

Conversely, consider an eigenfunction T/^pg) of -^2 with eigenvalue E]^2)- With equation 
0, we get H2ij'(2) = ^r^^^r2) = E]^2)i^h- Multiplying by At, we have, A^ Ai (A+V-^^s)) = 
E("2)^IV'"2)- With equation |10 



or from equation [10 

AtAi^l^ = Eli^^ly 
Applying Ai to equation |25|, AiA^ (^A^f-i/'fi)) = -^("i)^r^(i) with equation |l5 

H2 (Ar^[i)) = 



Hi = El) {Ati^l)) . 



(27) 



Equations ^ and imply that the hamiltonians Hi and H2 have identical eigenvalues 
except for the ground state -0°^-) of Hi (since Aiip'^^s^ = in equation ^ and this is 
unnormalizable) . In addition, we can see that if you know an eigenfunction of Hi, i.e. 
'0"]^), then an eigenfunction y4j~'0"^-) of H2 can be formed. Similarly, an eigenfunction Af'i/j^2) 
of Hi can be formed given an eigenfunction '0^^ of H2. The preceding analysis can then 
be extended to the hierarchy of hamiltonians discussed earlier. These observations are 
illustrated in figure |T] |T3|, |TB|. 

Herein lies a very important consequence of SUSYQM. The energy eigenfunctions of 
the hierarchy of hamiltonians are related by the SUSY operators A^. If one knows the 
eigenvalues and eigenfunctions of a particular Hn, then one can get the eigenvalues and 
eigenfunctions of its SUSY partner. 

Note that what we have discussed is SUSYQM in one dimension. There had been work 



done in doing SUSYQM in two or more dimensions 0, |I8|. To be able to apply one 
dimensional SUSYQM to the Hartmann potential, we will do a separation of variables of 
the resulting Schrodinger equation. To be able to form a hamiltonian from a separated one 
dimensional differential equation that can be an element of an Hgs, this one dimensional 
differential equation must be of the form of equation ^ (no first derivative term) and 
must yield an infinite tower of states as for Hn of figure |l] [|19|. As we will see, upon 
separation of variables for the Hartmann potential, only the radial equation yields an 
interesting SUSY. 
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3 Super symmetry in the Hartmann Hamiltonian 

With the SUSYQM concepts introduced in section we are now ready to obtain the 
eigenvalues and radial eigenfunctions of the Hartmann potential in spherical coordinates. 

The Schrodinger equation in spherical coordinates for a particle of mass fi subjected 
to the Hartmann potential in equation |1| is given by 



Assuming a solution 



i/j = i?(r)e(^)$(0) 
equation can be separated into three differential equations 

1 d^<l> 



-m 



d 



sinO dO 



+ 1) e = 



1 d 

I r' 

^2 



,dR' 
dr 



E + 



R = 



where 



2 I 2 2 

m + 7] a 



(28) 
(29) 

(30) 

(31) 
(32) 

(33) 



Looking at equations ^ to 0, we realize that these closely resemble the separated 
equations of the hydrogen atom |]2D|. As shown by reference [rU|, the only interesting 
separable SUSY in the hydrogen atom in spherical coordinates results from the radial 



equation. Their argument is as follows. Looking at equation 30, and comparing it with 



equation |T2|, we see that Vi = 0. Equation ^ on the other hand can be cast to a form 
similar to equation IT2| by multiplying it by a modulation factor /(™^^) ^ ^]^q eigenvalue 



[1— cos 

of f{cos9) is zero and no infinite tower of states can be generated. Hence, the $ and G 
solutions cannot be given by SUSYQM. They can be solved by conventional means 

1 



e(^) 



$(0) 



2n 
cos6), 



m 



0, ±1, ±2,... 



L = u'+\M\, 



0, 1, 2,... 



(34) 
(35) 



where (cosO) reduces to the associated Legendre polynomials when rja — >■ 0. 
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The radial equation |3^ can be cast into a form similar to O by letting 



r 



(36) 



Substituting equation ^ into equation yields 

1 £ L{L + 1) 7' 



2 dr-" 



+ 



with 



7 = 



2r2 



2 2 



u 



(37) 



(3^ 



Equation ^ is similar to that of the hydrogen atom's radial equation. We thus claim 
that we can obtain the eigenvalues and radial eigenfunctions by looking at the hamiltonian 
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n 1 L{L + 1) 7 1 / 7 

L 2dr^ 2r^ r 2\L + 1 



which yields a Ricatti equation (from equations |Tl| and |1^) 



L(L + 1) _7^ 1 



2r= 



7 



whose solution is 



r 2 VL + 1 
L + 1 



wl- 



dWL 
dr 



Equation ^ and ^ yield 



1 



+ 



7 



L + 1 



^ V2 V «r r L + 1 



(39) 



(40) 



(41) 



(42) 



tion 39 



From equation ^ and |T5|, we construct the SUSY-partner hamiltonian oiTii in equa- 

(43) 



n 



L+l 



7 



^ 2rfr2 2r=^ r 2 VL + 1 



Comparing equations and ^ and with equation we realize that 



7 



n 



L+l 



L+l 



2 VL + 1 

7 



(44) 



(45) 



2 VL + L 

Let us now start to build up the radial eigenfunctions and in the process get the 
eigenvalues. Given an \M\ value, the lowest L value is L = |M|, as can be seen in 
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equation It is apparent from equations ^ and ^ that we can build the states of the 
hierarchy of hamiltonians as in figure |I|. This is illustrated in figure 0. 

Since 'H\m\ {'H\m\+i) and Hl (i/i+i) differ only by a constant, (see equations |^ and 
^ ) every eigenfunction of 'H\m\ 0^\m\+i) will be an eigenfunction of (Hl^i). Hence, 
all we have to do is to solve for the eigenfunctions of Hl. The actual energy for Hl can 
be found by letting Hl act on the eigenfunctions of Hl- 

For an arbitrary L, equations and ^ give, for the ground states of TCl, '^(l)) the 
first order differential equation 



1 fd 
V2 \dr 



L + 1 



+ 



7 



L + 1 



which can easily be solved as 



where 



7 



L + 1 



(46) 



(47) 



(4^ 



Since L is arbitrary here, we realize that equation ^ is the expression for the eigen- 
function for the lowest rung (i.e. the ground state) of the tower of states for each of the 
hamiltonians in figure Since they are also eigenfunctions of Hl, we can write 



ul = Ml r^'^^exp{—KLr 



(49) 



where we illustrate them in figure I 



To get the actual energy, we let Hl of equation act on equation ^ 

1 d' 



HlUl 



+ 



2dr^ 

After some simplification, this yields 

El-- 



L{L + 1) 
2r2 



7 
r 



Ul 



fiEL 



A 



A 



2 4 I I 

V |eo| 



(50) 



(51) 



(^ + 1)- 

We characterize the energy by the L quantum number for the moment. From equation 
51|, we can label the energy levels of figure |^ as in figure ^. 

It is apparent from figure ^ that we have to label the solutions u as u\m\+i,\m\] 
1*1^^1+2,1^/1+1; "^lAfi+s, |M|+2; •••due to the energy of the states^ Knowing the eigenstates 
at the lowest rung of the hierarchy of hamiltonians, M|Af|+i, |Af|; "^11/1+2,11/1+1; W|a/|+3, |a/|+2; 
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etc., we can determine the other states by the action of on these eigenstates as in 
figure |I|. This is illustrated in figure |^. 

Note that, for instance, u\m\+3,\m\] ^i|M|+3,|M|+i; ^i|M|+3,|M|+2; • • - have the same energy 
~(Tm1+3F' similarly for other states at the same energy level. It is then evident that 
given > \M\ + 1, un,n-i', ujsi^n-i] un,n-z] ■ ■ ■ ] un,\m\ will all have the same energy 
— Hence, we can say that 



E 



N 



-A, A = r/V|eo|, iV = L + l + n', n' = 0,l,2,... (52) 



which means that the energy is actually labelled by and not L. Equation ^ agrees 
with reference [Q. 

From figure |^, equations ^ ^ ^ and it can be shown that 



M\+l,\M\\f) 



-R|Af|+2,|Af|+l('^) 
-R|M|+2,|M|(''') 




X 



|Af|+3/2 



|Af|+5/2 



I A/ 1 +3/2 



|Af|+2j 
2|M| +2 



r(2|M| +3) 

1 

r(2|M| +5) 
1 



1/2 



^|A/|g-7r/[|Af|+l] 



1/2 



^|Af|+lg-7r/[lAf|+2] 



2(|Af|+2)r{2|Af|+3) 
2'yr 



1/2 



^lAf|g-7r/[|A/l+2] 



|Af|+2 



(53) 

(54) 
(55) 



etc. where Rnl are normalized by 



\R 



NL 



1^ r'^dr 



and the gamma function properties used are [21 



roo 

T{z) = / e-H'-^dt; T{z + 1) = zT{z) 
Jo 



(56) 



(57) 



Equations ^ ^ and ^ agree with the normalized Rnl of reference 0. It is obvious 
that we can eventually get all the expressions of Rnl- 

To illustrate the procedure more concretely, let us show how we can get -R|A/|+2,|Af|('^)- 
From figure ^, 



M|Af|+2,|Af| ~ ^l^vfl «|Af|+2,lM|+l-' 



(5^ 



^ Note that in equation 2 and in some of the foUowing equations, we use ~ instead of = since the 
procedure outhned here cannot automaticaUy normahze the eigenf unctions. At the end, we normahze 
the eigenfunctions using equation pq. 
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Since «|M|+2,jA/|+i is at the lowest rung of H\m\+i, we can use equations ^ and ^ which 
give 



u\M\+2,\M\+i ^ r' e ' 



From equations ^ 2 and we get 

f d \M\ + 1 



U\M\+2,\M\ 



7 



dr 



1^1 + 1, 



^|M|+2g-7r/[|M|+2] 



which leads to 



W|M|+2,|M| ~ re 



-yr/[\M\+2] _JM\ 



7r 



|A/|+1 



(|M|+2)(|M| + 1) 



Using equations 61 and 36 and rearranging terms, we get 



i?|Mi+2,iMi ~ -e-^^/[l*^l+^l rl^-^l 2 + 2|M| 



2^r 



\M\ + 2 



or 



Rm.2m = -ATe--/!!-!-^] rl-l (2 + 2 |M| - 



(59) 



(60) 



(61) 



(62) 



(63) 



where A/" is the normalization constant. Normalizing equation |6^ using equations ^ and 



57 leads to equation 55 
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4 Conclusion 



The ring-shaped Hartmann potential was first introduced in 1972 ||T|. It has been 
used to describe ring-shaped molecules like the benzene molecule in theoretical chemistry 
p[. In an attempt to introduce the methods of supersymmetric quantum mechanics in 
quantum chemistry, we obtain the eigenvalues and radial eigenfunctions of the Hartmann 
potential in spherical coordinates using SUSYQM techniques. 

A key result in SUSYQM is the intimate relationship of the eigenvalues and eigenfunc- 
tions of the hierarchy of SUSY-partner hamiltonians. This can be very useful in solving 
the Schrodinger equation of a complicated hamiltonian if its SUSY-partner hamiltonian 
is easily solvable. 

A very useful result in the present discussion is the fact that A~ifj^ = 0. This enabled 
us to solve a first order differential equation (as in equation ^Gj) instead of the second 
order Schrodinger differential equation to give us the eigenfunctions and eigenvalues of 
the states at the lowest rung of the tower of states of each of the hamiltonians in the 
hierarchy. The rest of the eigenfunctions and eigenvalues are then easily computed by 
applying the corresponding Aj^ operators to these eigenfunctions. 

As already indicated in this paper, the separated equations of the Hartmann potential 
and the hydrogen atom greatly resemble each other. In fact, the technique employed 
here was first applied to the hydrogen atom [|T3|. A number of studies of the SUSY 
features of the Coulomb problem in hydrogenic atoms have been made over the past years 
[p!9| , 0, |2^, ^ . These studies may very well lead to some further insights into the workings 
of SUSY in the Hartmann potential due to the similarity of its separated equations with 
that of the hydrogen atom. 

The preceding paragraph actually led the author to believe that if the Hartmann 



problem is treated using the SUSY formulation of Haymaker and Rau p4[ , the eigenstates 
of different values of ija"^ can be related. This will be developed in a future publication. 

With the above considerations, the author believes that quantum chemistry is a field 
where the concepts and techniques of SUSYQM can be put to good use. 
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Figure Captions 



1. The hierarchy of hamiltonians and the action of the operators on the degenerate 
eigenstates 

2. The energy states of the hierarchy of SUSY-partner hamiltonians from the Hartmann 
potential 

3. The hierarchy of hamiltonians of the Hartmann potential and their ground states. 
The Hl here are the actual radial hamiltonian for a particular L value. 

4. Figure 3 with the energy levels labelled. 

5. The energy eigenstates of the Hartmann potential. The action of the A'l operators 
are explicitly shown to indicate how the other states are obtained from the states 
at the lowest rung of the hierarchy of hamiltonians. 
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Figure 1: The hierarchy of hamiltonians and the action of the operators on the 
degenerate eigenstates 
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7i|M| Ti-lMl+l 'H\M\+2 — 

Figure 2: The energy states of the hierarchy of SUSY-partner hamiltonians from the 
Hartmann potential 
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U\M\+2 



U\M\+1 



U\M\ 

H\M\ H\M\+1 H\M\+2 

Figure 3: The hierarchy of hamiltonians of the Hartmann potential and their ground 
states. The Hi^ here are the actual radial hamiltonian for a particular L value. 



17 



H\M\ H\M\+1 H\M\+2 



Figure 4: Figure 3 with the energy levels labelled. 
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(|M|+3) 



^|M|+3,|M| 



^\M\ ^\M\+l 
1i|M|+3,|M|+l M|M|+3,|M|+2 



-A U\M\+2,\M\ 

(lM|+2)^ 



4 + 
^|M| 



'»|M|+2,|M|+1 



-A M|M|+1,|M| 



H\M\ H\M\+l H\M\+2 



Figure 5: The energy cigcnstatcs of the Hartmann potential. The action of the A]^ 
operators are exphcitly shown to indicate how the other states are obtained from the 
states at the lowest rung of the hierarchy of hamiltonians. 
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